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ABSTRACT  In  a  horizontal  layer  of  fluid  heated  below  and  cooled 
above,  cellular  convection  with  horizontal  length  scale  comparable 
to  the  layer  depth  occurs  for  small  enough  values  of  the  Rayleigh 
number.  As  the  Rayleigh  number  is  increased,  cellular  flow 
disappears  and  is  replaced  by  a  random  array  of  transient  plumes. 
Upon  further  increase,  these  plumes  drift  in  one  direction  near  the 
bottom,  in  the  opposite  direction  near  the  top  of  the  layer  with 
the  axes  of  plumes  tilted  in  such  a  way  that  horizontal  momentum  is 
transported  upwards  via  the  Reynolds  stress.  With  the  onset  of 
this  large  scale  flow,  the  largest  scale  of  motion  has  increased 
from  that  comparable  to  the  layer  depth  to  a  scale  comparable  to 
the  layer  width.  The  conditions  for  occurrence  and  determination 
of  the  direction  of  this  large  scale  circulation  are  described. 
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INTRODUCTION  This  is  a  report  on  a  laboratory  study  of 
convection  in  a  horizontal  layer  of  fluid  uniformly  and  steadily 
heated  from  below  and  cooled  from  above.  This  is  the  arrangement 
which  leads  to  cellular  convection  for  a  range  of  values  of  the 
Rayleigh  number  (the  dimensionless  measure  of  temperature 
difference,  defined  below),  sufficiently  small  but  above  a 
critical  value.  The  horizontal  scale  of  these  cells  is  comparable 
to  the  depth  d  of  the  layer.  At  successively  larger  values  of  this 
parameter,  a  number  of  transitions  in  the  flow  pattern  as  well  as 
in  the  heat  flux  are  observed  (1-6).  Most  of  these  changes  are 
within  the  regime  of  cellular  flows.  We  shall  describe  here  a 
further  transition  which  leads  to  a  very  different  scale  of  motion 
and  very  different  transport  properties.  In  this  regime  the  flow 
is  no  longer  cellular.  It  is  a  flow  having  primarily  two  scales  of 
motion.  The  smaller  scale  flow,  with  horizontal  scale  comparable 
to  the  layer  depth  d,  is  best  described  as  transient  bubbles  or 
plumes  that  have  an  organized  tilt  away  from  the  vertical.  The 
larger  scale  flow,  having  scale  l,  is  a  horizontal  flow  with 
vertical  shear  such  that  the  flow  is  oppositely  directed  near  the 
bottom  and  the  top.  L  is  usually  the  container  width  which  in 
these  experiments  is  an  order  of  magnitude  larger  than  d.  The 
large  scale  flow  is  apparently  maintained  against  viscous 
dissipation  by  the  Reynolds  stress  divergence  of  the  tilted  plume 
motions. 
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APPARATUS  &  PROCEDURE  The  convecting  fluid  occupied  a  space  of 
depth  d  and  horizontal  extent  L  by  L,  where  d/L  was  a  small  number 
of  order  10~2  to  10~1.  Most  of  the  new  results  reported  below  were 
obtained  with  d  =  2  cm  and  5  cm,  L  =  48  cm.  The  fluid  layer  was 
bounded  above  and  below  by  metal  blocks  whose  the'rmal  conductivity 
was  several  thousand  times  that  of  the  liquids  used.  Even  in  the 
most  extreme  case  within  this  set  of  experiments,  when  the 
convective  heat  flux  was  thirty  times  the  conductive  heat  flux,  the 
boundaries  had  conductivity  two  orders  of  magnitude  higher  than  the 
effective  conductivity  of  the  fluid  layer.  Thus  the  boundary 
condition  is  one  of  constant  temperature  along  both  top  and  bottom 
boundary.  On  all  its  lateral  boundaries  the  fluid  layer  was  bounded 
and  insulated  by  3  cm  of  Plexiglas  and  5  cm  of  styrofoam.  Normally 
the  apparatus  was  levelled  so  that  its  slope  was  no  more  than  .0006 
cm  in  30  cm.  The  flatness  and  parallelism  of  the  horizontal 
boundaries  were  comparable  to  this.  The  fluids  used  were  water  and 
some  of  the  silicone  oils. 

The  auxiliary  equipment  for  control  and  measurement  of  heat 
flux  are  described  briefly  here,  and  in  detail  in  (3).  In  addition 
to  the  convecting  layer,  six  layers  of  solid  materials,  each  51  cm 
square,  were  used  in  a  stack.  The  bottom  layer  was  a  10  cm  thick 
aluminum  block  with  the  heater  (made  of  a  fine  mesh  of  resistance 
material)  attached  to  its  bottom.  Above  this  lay  a  0.64  cm  thick 
Plexiglas  layer  and  above  this  a  5  cm  thick  copper  layer.  .Above 
this  lay  the  convecting  fluid.  Above  the  fluid  were  5  cm  of 
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copper,  then  0.64  cm  of  Plexiglas,  then  a  10  cm  aluminum  block 
containing  channels  for  circulating  cooling  fluid.  (In  some  of  the 
earlier  experiments,  2.5cm  thick  aluminum  blocks  were  used  in  place 
of  the  copper  blocks.)  The  low  conductivity  Plexiglas  gives  rise  to 
measurably  large  temperature  differences  between  the  metal  blocks 
and  was  used  to  measure  the  total  heat  flux,  (a  method  devised  by 
Malkus  (l)i.  The  entire  stack  was  placed  in  a  Plexiglas  tank, 
which  was  surrounded  on  all  sides  by  styrofoam  of  5  cm  thickness. 

The  procedure  mi  to  set  the  heater  voltage,  then  adjust  the 
temperature  of  the  cooling  water  so  that  when  the  temperatures 
measured  in  the  four  metal  blocks  reached  steady  state,  the  mean 
temperature  of  the  converting  fluid  remained  at  room  temperature. 
Thus  the  total  heat  flux  H  was  prescibed  (although  not  measured)  by 
the  input  voltage.  The  temperature  is  a  constant  on  the  top  and 
bottom  boundaries  of  the  fluid,  but  the  magnitude  of  the  difference 
AT  between  these  constants  is  determined  indirectly  through  the 
convective  heat  flux  of  the  fluid. 

The  Rayleigh  number  was  computed  from  the  measured  AT.  The 
total  heat  flux  H  through  the  fluid  wa s  equated  to  the  conductive 
flux  through  either  of  the  Plexiglas  layers.  The  conductivity  of 
the  Plexiglas  was  determined  in  terms  of  the  known  conductivity  of 
water,  by  performing  the  experiment  with  a  0.5  cm  depth  of  water 
and  AT  smaller  than  the  critical  value  for  onset  of  convection. 

To  study  the  geometry  of  the  time-dependent  flow,  the  following 
arrangement  was  used  to  obtain  (x,t)  photographs  from  which 
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information  about  the  Eulerian  velocity  could  be  deduced.  Here  x 
refers  to  the  horizontal  coordinate  along  a  fixed  line  within  the 
fluid,  perpendicular  to  the  optic  axis  of  the  camera,  and  t  refers 
to  the  time.  We  also  let  y  be  the  coordinate  along  the  optic  axis 
and  z  the  vertical  coordinate.  The  fluid  contains  suspended  tracer 
particles  which  are  small  light  reflecting  flakes  that  become 
aligned  by  the  shear  of  the  flow.  The  illuminated  region  is  2  to  3 
mm  in  diameter,  48  cm  long  along  x,  and  was  chosen  to  lie  either 
near  the  bottom  or  the  top  of  the  fluid  layer.  The  lights  used 
were  either  two  2-watt  zirconium  arc  lamps  or  a  5  milliwatt 
continuous  laser.  With  the  illuminated  region  fixed  in  space,  a 
wedge  was  moved  under  the  camera  in  such  a  way  that  the  camera 
would  rotate  about  a  horizontal  axis  (parallel  to  the  x  axis) 
through  its  lens.  The  resulting  photograph  shows  x  along  the 
abscissa,  t  along  the  ordinate.  For  example,  if  x  were  a  line  near 
the  bottom  of  the  layer,  perpendicular  to  the  longitudinal  axis  of 
roll  convection  cells,  the  (x,t)  photograph  would  be  as  shown  in 
figure  2a.  A  tracer  particle  which  is  within  the  illuminated 
region  and  properly  oriented  so  as  to  send  light  into  the  camera 
produces  a  mark  on  the  film.  Individual  particles  generally  remain 
so  for  only  a  small  fraction  of  the  total  exposure,  but  long  enough 
to  make  a  short  line  segment  in  the  picture,  as  is  clearly  seen  in 
figure  2a.  The  local  slope  dx/dt  of  the  direction  field  in  the 
(x,t)  photograph  produced  by  these  line  segments  gives  the 
horizontal  velocity  component  as  a  function  of  x  and  t.  Cell 
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boundaries,  from  which  fluid  is  diverging  or  toward  which  it  is 
converging,  are  readily  identified,  not  only  from  this  divergence 
or  convergence  but  from  changes  in  the  general  brightness  or 
darkness  brought  about  by  the  tendency  of  the  plates  to  become 
oriented  in  regions  of  high  shear.  In  steady  convection  as  in  2a, 
cell  boundaries  remain  fixed  for  all  time. 

Stray  heat  sources  that  might  produce  lateral  temperature 
gradients  were  minimized  by  the  physical  arrangement  of  the 
equipment.  The  convection  tank  along  with  the  photographic 
apparatus  was  enclosed  in  a  light-tight  house.  Except  for  the 
lights,  the  motor  drive  for  the  camera  (and  sometimes  the 
experimental ist)  and  of  course  the  heater  for  the  fluid  layer,  all 
other  heat  producing  instruments  of  control  or  measurement  were 
kept  outside  the  house. 

After  the  large  scale  flow  was  detected,  several  additions  to 
the  apparatus  were  made  to  test  what  controlled  the  direction  of 
this  circulation.  The  first  was  a  paddle  made  of  0.32  cm  diameter 
stainless  steel  rod,  47  cm  long,  which  could  be  drawn  in  the 
x-direction  perpendicular  to  its  length.  The  second  addition  was  a 
jack  to  tilt  the  entire  apparatus.  The  third  was  the  addition  of  a 
heater  in  the  side  boundary  of  the  fluid  layer.  A  35  cm  length  of 
nichrome  wire  was  embedded  in  a  0.64  cm  thick  sheet  of  Plexiglas. 
This  sheet  replaced  a  similar  one  which  was  part  of  the  side  wall. 
The  heater  wire  occupied  a  horizontal  position  approximately  one  cm 
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higher  than  the  bottom  of  the  convecting  layer  and  1.3  cm  outside 
the  fluid. 

Finally,  the  large  scale  circulation  was  also  observed  in  a 
cylindrical  annular  tank  of  fluid,  uniformly  heated  below  and 
cooled  above.  The  diameter  of  the  outer  Plexiglas  cylinder  was 
45.8  cm,  the  gap  width  between  the  two  cylinders  was  4.3  cm,  the 
fluid  depth  d  was  10.6  cm.  The  lid  was  floated  on  an  air  bearing, 
but  restrained  from  rotating  except  when  testing  the  effect  of 
imposed  motion  of  the  lid  upon  direction  of  large  scale 
circulation.  For  cooling  the  top  boundary,  the  lid  carried  a 
volatile  fluid  open  to  the  air. 

OBSERVATIONS  Figure  1  summarizes  some  observations  of  the  flow 
found  in  the  various  ranges  of  the  parameters.  These  parameters 
are  the  Rayleigh  number  R  and  the  Prandtl  number  P  defined  as 
fol 1 ows: 

R  =  9*  ATd3  ,  P  =  ^ 

<V  K 

where  g  is  the  acceleration  of  gravity,  a  the  thermal  expansion 
coefficient,  <  the  thermal  diffusivity,  v  the  kinematic  viscosity, 
d  the  depth  of  the  fluid  layer,  and  AT  the  temperature  difference 
between  bottom  and  top  boundaries.  The  new  observations  refer  to 
flows  found  to  the  upper  left  of  the  curve  labelled  V  in  figure  1. 
To  the  lower  right  of  .his  curve  the  flows  are  cellular:  as 
labelled,  the  flow  is  steady,  two-dimensional  in  the  parameter 
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range  between  curves  I  and  II.  It  is  steady  three-dimensional 
cellular  between  curves  II  and  III,  and  time-dependent 
(quasi-periodic) ,  three  dimensional,  and  cellular  between  curves 
III  and  V.  8y  cellular  we  mean  that  a  fluid  parcel  initially  at 
some  horizontal  location  is  always  confined  to  the  vicinity  of  that 
location,  horizontal  excursions  being  limited  to  lie  between  cell 
boundaries  which  have  horizontal  separation  comparable  to  the  layer 
depth.  (This  is  certainly  so  between  curves  I  and  III  and  at  least 
largely  so  between  III  and  V.)  Since  most  of  these  states  are 
history  dependent,  the  actual  transition  point  varies  depending 
upon  how  it  was  attained.  Figure  1  obtains  when  the  state  of  the 
system  was  reached  from  lower  values  of  Rayleigh  number,  by  small 
increments  in  R.  Further  explanation  can  be  found  in  (5). 

The  curve  V  was  defined  in  (2)  as  the  Rayleigh  number  above 
which  irregular  fluctuations  in  temperature  at  a  point  mid-depth  in 
the  fluid  layer  were  nearly  always  present.  Although  somewhat 
subjective,  it  is  a  useful  criterion. 

One  of  the  major  differences  between  the  flows  occurring  to  the 
upper  left  of  V  from  that  to  the  lower  right  is  seen  in  the  (x,t) 
photographs  of  figure  2.  Both  2b  and  2c  show  a  line  in  the  fluid 
layer  0.2  cm  above  the  bottom  boundary.  Figure  2b,  at  R  =  10®, 

P  =  .86  X  10^  is  for  parameter  values  to  the  right  of  curve  V. 

Figure  2c,  at  the  same  Rayleigh  number,  =  10®,  but  P  =  7,  is  to 
the  left  of  7.  In  2b  there  are  cell  boundaries  which,  while  they 
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may  oscillate  laterally,  are  always  identifiable.  In  2c  there  are 
no  cell  boundaries.  The  former  is  a  cellular  flow  with  bright  or 
dark  regions  moving  repeatedly  from  one  cell  boundary  to  a 
neighboring  one.  These  regions  have  been  shown  in  (5)  to 
correspond  to  hot  or  cold  spots  convected  around  the  cell.  In  2c, 
hot  bubbles  or  plumes  off  the  bottom,  randomly  spaced  in  position 
and  time,  form  and  vanish.  Cold  bubbles  off  the  top  likewise  form 
and  vanish,  much  as  described  by  the  bubble  model  (3).  At  any 
instant  of  time,  their  horizontal  soacing  is  approximately  d  or 
somewhat  less.  It  should  be  pointed  out  also  that  figure  2c  is  not 
an  initial  transient  state;  it  is  the  statistically  steady  state 
that  is  observed  one  hour  or  100  hours  after  starting  convection. 

At  P  =  7,  R  =  5  X  104  the  flow  is  similar  to  2b,  so  that 
increasing  R  at  fixed  P  also  leads  to  the  change  of  flow  depicted 
when  going  from  figure  2b  to  2c. 

The  other  major  difference  -the  one  to  which  we  particularly 
wish  to  call  attention  in  this  report  -is  seen  in  figure  2d  which 
shows  an  (x,t)  photograph  for  P  =  7,  R  =  2.4  X  10®,  taken  0.2  cm 
above  the  bottom  boundary.  It  shows  that  the  plumes  have  a  net 
drift  from  left  to  right  with  an  average  speed  of  7  X  10~2  c;n 
sec"1.  A  similar  photograph  taken  0.2  cm  below  the  top  showed 
plumes  drifting  from  right  to  left  with  the  same  average  speed. 

From  the  (x,t)  photos  it  is  not  known  if  this  apparent  motion  is 
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a  wave  propagation,  or  if  there  is  an  associated  mass  flux.  Thymol 
blue  was  added  to  the  water  so  that  neutrally  buoyant  dye  could  be 
produced  at  an  electrode,  and  the  dyed  fluid  was  followed  with 
time.  Figure  3  shows  a  sequence  at  30  sec  intervals,  of  a  vertical 
slice  (x,z)  view  of  the  fluid.  In  this  case  the  layer  depth  was  10 
cm.  The  dye  front  along  the  bottom  progressed  at  an  average  speed 
u  of  3  X  10-2  cm  sec"*  from  left  to  right,  and  actually  reached 
the  opposite  wall,  by  which  time  it  was  quite  diffuse.  These 
magnitudes  are  not  easily  compared  from  case  to  case  since  only  the 
x-component,  u,  is  measured,  whereas  in  some  repetitions  of  the 
experiment  there  was  a  velocity  component  v  in  the  y-di recti  on  as 
well.  Maximum  values  of  TT  on  the  order  of  10-1  cm  sec-1  have  been 
observed.  In  figure  3,  as  the  dyed  fluid  becomes  buoyant,  it  rises 
into  the  interior  of  the  layer  where  it  is  seen  to  become  entangled 
in  the  motion  of  the  plume.  It  often  rises  to  become  caught  in  the 
leftward  drift  at  the  top  of  the  layer.  In  figure  3,  a  general 
tilt  of  the  plumes  from  lower  right  to  upper  left  can  be  discerned 
even  though  the  flow  is  quite  turbulent.  Thus  we  conclude  that 
there  is  not  only  a  (time-and/or  horizontal)  mean  Eulerian 
velocity  Tf(z) ,  but  a  net  horizontal  Lagrangian  transport, 
extending  the  entire  width  of  the  tank.  Observations  of  the 
occurrence  of  cells,  plumes,  and  large  scale  flow  in  various 
regions  of  parameter  space  are  summarized  in  Figure  1. 
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On  repetitions  of  the  experiment  the  large  scale  flow  was  seen 
to  be  sometimes  from  left  to  right  along  x  at  the  bottom  and  from 
right  to  left  along  x  at  the  top  (U  counterclockwise).  On  other 
realizations  the  flow  was  clockwise.  On  still  others,  the  large 
scale  flow  was  in  the  y-direction  (7)  and  either  clockwise  or 
counterclockwise.  It  may  also  be  at  some  other  angle  to  the 
x-direction.  When  the  flow  was  "u  and  clockwise,  this  flow  was 
often  nearly  independent  of  y,  but  on  another  repetition  ic  ; ght 
be  that  the  center  of  the  large  scale  rising  was  at  the  left  wall 
at  some  value  of  y,  but  at  other  values  of  y,  the  center  of  the 
rising  wandered  to  other  values  of  x. 

When  this  large  scale  flow  was  first  discovered,  various  tests 
were  made  to  see  if  it  were  inadvertently  being  driven  directly  by 
horizontal  asymmetries.  The  levelling  of  the  apparatus  was 
checked.  Horizontal  temperature  gradients  in  the  boundaries  were 
measured  to  check  the  possibility  of  malfunction  of  the  cooling  or 
heating  system.  In  the  aluminum  blocks  with  the  heater,  or  the 
cooler,  horizontal  gradients  were  of  order  10“^  °C  across  35  cm. 

In  the  middle  two  (copper)  blocks,  which  do  not  contain  any  thermal 
forcing,  the  horizontal  gradients  were  larger:  10"?  °C  in  25  cm. 
However,  this  magnitude  of  horizontal  gradient  would  drive,  in 
water,  a  horizontal  velocity  u  =10"^  cm  sec'^-  which  is  two  orders 
of  magnitude  smaller  than  observed.  Although  horizontal 
asymmetries  must  certainly  be  present,  due  to  errors  in  levelling 
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the  apparatus,  or  due  to  small  differences  in  temperature  of  the 
Plexiglas  side  walls,  these  had  no  observable  effect  at  lower 
Rayleigh  number,  or  even  at  the  same  Rayleigh  number  (with  the  same 
ratio  of  side  wall  area  to  horizontal  area),  when  a  higher  Prandtl 
number  fluid  was  used  and  no  large  scale  flow  was  observed. 

Convinced  that  the  large  scale  flow  was  not  externally  driven, 
we  attempted  to  check,  as  to  order  of  magnitude,  that  the  large 
scale  flow  and  the  tilt  of  the  plumes  were  consistent  with  the 
Reynolds  equations  relating  a  (time-  and  horizontal-)  mean  flow 
with  the  Reynolds  stress  associated  with  fluctuations.  The 
principal  equation  here  is  the  mean  horizontal  momentum  equation, 
which  if  the  tank  were  horizontally  unbounded  (or  in  the  case  of 
the  cyl indrical  annulus)  is: 

(uV)z  =  vTTzz  ;  (*) 

the  overbar  indicates  a  horizontal  and  time  average,  u'  and  w'  are 
the  fluctuations  from  the  mean  and  v  is  the  kinematic  viscosity. 
(With  the  finite  horizontal  extent  there  could  conceivably  also  be 
a  non-zero  mean  horizontal  pressure  gradient  term  in  this  equation 
--  which  indeed  might  well  be  important  were  there  external  driving 
--  but  we  omit  this  here.)  Integrating  vertically  from  the  bottom 
wall  we  obtain  u 1  w’ 1  =  v(ITz  -  "uz),  "uz  being  the  mean  shear  at  the 
wall.  Now  for  a  rough  estimate  we  may  perhaps  suppose  that  the 
mean  flow"u  has  a  profile  in  z  something  like  u=-U  sin(2irz/d), 
taking  the  origin  of  z  at  the  center  of  the  fluid  layer  so 
-d/2<z<d/2.  This  would  correspond  to  a  Reynolds  stress 
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u1  w'  =  (4itUv/d)  cos2(irz/d).  Thus  with  U>0,  corresponding  to  flow 
from  left  to  right  in  the  lower  half  of  the  layer  (as  in  fig.  3), 
u‘ w1  must  be  negative  and  its  maximum  magnitude,  achieved  at  the 
center,  should  be  about  10v/d  times  the  maximum  of  U.  Detailed 
measurements  of  the  fluctuating  velocity  components  u1  and  w',  and 
of  their  correlation,  have  not  yet  been  made.  However  visual 
observation  and  photographs  like  fig.  3,  together  with  the  (x,t) 
photographs,  indicate  that 

a)  the  fluctuations  are  on  the  whole  associated  with  motions  in 
the  tilting  plumes; 

b)  the  latter  tilt  in  such  a  way  (about  45°  from  lower  right  to 
upper  left  in  cases  like  figure  3)  that  u  V  has  opposite  sign  to 
that  of  H  in  the  lower  half,  and  w'  is  characteristically  of 
about  the  same  magnitude  as  u * ; 

c)  the  maximum  magnitudes  of  u',  w' ,  and  U  are  all  about  0.1 
cm/sec  in  our  experiments  in  which  the  large  scale  flow  is 
observed. 

We  may  thus  estimate  u'w'  as  about  Q.2X10"2  cm^sec"^  if  we 
include  factors  of  about  1/2  for  the  x-average  (assuming  u‘  and  w' 
have  sinusoidal  variation  with  x)  and  also  for  the  time  average. 

The  last  is  a  very  rough  estimate  from  the  fraction  of  the  time 
that  plumes  are  seen  at  any  location  on  the  (x,t)  photographs. 

With  d=5  cm  and  v=.01  cm^sec"^  we  get  (10v/d)umax5  0.2X10*2  cm2 
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sec‘2,  so  there  does  not  appear  to  be  any  clear  inconsistency  with 
the  Reynolds  equation  (*).  One  may  regard  the  mean  flow  as  being 
"driven"  by  the  Reynolds  stress  associated  with  the  tilting  plumes. 
We  emphasize  that  such  a  statement  is  a  description,  but  not  an 
explanation,  of  the  large  scale  flow.  The  Reynolds  equation  would 
also  be  consistent  with  no  mean  flow,  and  no  tilting  of  the  plumes 
--  and  in  fact  this  appears  to  be  what  occurs  below  a  Rayleigh 
number  of  about  2X10^  in  water.  The  theoretical  understanding  of 
why  this  transition  occurs,  and  why  at  R  =  2X10^,  remains  a 
challenge.  The  phenomenon  seems  to  resemble  the  onset  of  steady 
convection  in  that  the  physically  realized  "state"  changes  to  one 
of  a  set,  each  member  of  which  individually  has  less  symmetry  than 
the  equations  and  boundary  conditions  describing  the  motion.  Such 
symmetry-breaking  occurs  in  many  mathematical  models  involving 
changes  of  stability  of  stationary  or  periodic  solutions  and 
associated  bifurcations,  and  Rayleigh's  explanation  of  the  onset  of 
cellular  convection  in  such  terms  has  been  followed  in  spirit  by 
essentially  all  of  the  many  subsequent  explorations  of  other 
aspects  of  cellular  convection.  Whether  such  concepts  as  stability 
and  bifurcation  can  usefully  contribute  to  an  understanding  of  a 
transition  entirely  in  the  turbulent  range  seems  less  clear;  this 
report,  in  any  case,  is  about  the  observation  of  the  transition, 
not  its  theoretical  interpretation. 

What  determines  the  direction  of  the  large  scale  flow?  We  were 
not  successful  in  altering  the  direction  of  the  established  large 
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scale  flow  by  drawing  the  paddle  along  the  bottom  of  the  convecting 
layer,  although  we  tried  various  speeds  including  one  approximately 
the  speed  of  the  large  scale  flow.  We  did  not  attempt  this  while 
the  Rayleigh  number  was  being  increased  through  R  »  2  X  10®.  It 
was  not  possible  to  match  the  time  scales  involved.  Similarly,  in 
the  annular  tank  it  was  not  possible  to  alter  the  direction  of  the 
already  established  large  scale  flow,  by  rotating  the  lid  through 
one  revolution,  at  various  speeds,  in  either  the  clockwise  or 
counterclockwise  direction. 

When  the  entire  convection  tank  was  tilted  by  various  angles 
up  to  a  maximum  of  0.23  degrees,  neither  the  already  established 
large  scale  flow  nor  the  large  scale  flow  that  set  in  as  R  was 
increased  from  zero  to  beyond  2  X  10®,  was  determined  by  the  tilt 
of  the  layer.  Certainly  it  is  known  that  if  the  tilt  angle  is 
large  enough  the  direction  of  the  flow  is  determined  by  that  tilt, 
but  for  the  tilts  used  here  (which  were  10^  times  larger  than  any 
tilt  of  the  layer  due  to  levelling  error)  there  was  no  such 
determination. 

When  the  fluid  layer  was  heated  from  the  side,  the  direction  of 
the  large  scale  flow  could  be  determined  if  the  heating  rate  was 
large  enough.  Several  different  heating  rates  were  tried.  For 
example,  when  the  power  to  the  side  heater  was  P  =  3  watts,  this 
heating  rate  did  not  reverse  an  already  established  large  scale 
flow.  The  rising  motion  induced  at  the  heated  side  wall  was 
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accompanied  by  sinking  motion  a  short  distance  (less  than  d)  a way, 
and  the  rest  of  the  layer  behaved  as  without  the  side  heating. 

When  the  side  heater,  at  P  =  3  watts,  was  on  while  R  was  increased 
from  zero  through  2  X  10®,  the  direction  of  the  large  scale  flow 
was  not  determined  by  the  side  on  which  the  heater  was  placed. 
However,  when  the  heating  rate  was  P  =  11  watts,  as  R  was  increased 
from  zero  to  2  X  10®,  the  large  scale  flow  was  "u,  clockwise,  if  the 
heater  was  on  the  left,  and  a  "u,  counterclockwise,  if  the  heater 
was  on  the  right.  Furthermore,  with  P  =  11  watts,  an  established 
large  scale  flow  could  be  made  to  reverse  direction  by  switching 
the  heater  to  the  other  side  of  the  tank.  The  response  was  within 
an  hour;  the  thermal  diffusion  time  of  the  Plexiglas  spacer  is 
approximately  one-half  hour.  Also,  after  the  side  heater  was 
turned  off,  the  direction  of  the  large  scale  flow  remained 
unchanged  for  several  days  thereafter.  This  was  with  R  =  107 .  The 
heating  rate  of  11  watts  caused  an  increase  in  temperature  of 
approximately  1°C  at  the  outside  of  the  Plexiglas  spacer,  about  1 
cm  from  the  fluid.  The  temperature  difference  in  the  vertical  was 
2°C.  This  heating  alone,  with  the  vertical  Rayleigh  number  R  =  0, 
produced  a  large  scale  overturning  with  horizontal  velocity  u(x) 
measured  at  a  distance  z  =  0.2  cm  below  the  top  of  the  5  cm  deep 
layer.  The  measurement  was  conveniently  made  from  the  slopes  of 
particles  in  an  (x,t)  photograph.  The  maximum  value  of  u,  which 
occurs  at  x  =  14  cm  from  the  heated  wall,  for  z  =  0.2  cm,  was  8  X 
10"3  cm  sec"l.  This  is  still  one  order  of  magnitude  less  than  the 
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u  produced  by  the  tilting  plumes. 

Finally  large  scale  flow  was  observed  in  the  cylindrical 
annulus,  showing  that  a  pressure  difference  between  two  lateral 
boundaries  is  not  necessary  to  produce  this  flow.  The  small  scale 
plumes  are  considerably  narrower  than  the  depth  d  of  the  layer,  but 
they  tilt  approximately  45°  and  in  this  sense  they  occupy  a  hori¬ 
zontal  distance  of  about  d.  The  large  scale  flow  was  observed  to 
continue  around  the  annulus,  in  one  direction  near  the  bottom,  and 
in  the  opposite  direction  near  the  top.  With  the  geometry  thus 
unbounded  in  x,  the  large  scale  flow  has  horizontal  wave  number 
zero. 

Movies  of  the  flow  in  the  annulus  and  of  an  (x,z)  slice  in  the 
square  tank  are  available. 

The  heat  flux,  non-dimensional i zed  to  equal  the  Nusselt  number 
times  Rayleigh  number,  is  plotted  against  Rayleigh  number  in  figure 
4.  (The  Nusselt  number  is  the  ratio  of  total  to  conductive  heat 
flux.)  Discrete  changes  in  slope  were  seen  at  the  Rayleigh  numbers 
Rin  =  3.3  X  104,  Riv  =  5.5  X  104,  Ry  =  1-0  X  105,  Ryi  =  5.0  X  105, 

Ryn  =  1*1  x  10®,  the  ]ast  three  of  whlch  are  shown  1n  fig*  4> 

Some  of  these  are  in  good  agreement  with  the  Malkus  heat  flux 
transitions  Rjjj  =  3  X  104 ,  RIy  =  5.5  X  104,  Ry  =  1.7  X  10®,  Ryj  = 
4.1  X  105,  Ryu  =  8.5  X  105,  Ryn  =  1.6  X  IQ6.  It  appears  that  the 
onset  of  the  large  scale  circulation  might  be  associated  with  the 
slope  change  at  Ryu  =  1.1  X  10®. 


DISCUSSION  As  the  Rayleigh  number  is  increased,  cellular  flow 
vanishes  and  is  replaced  by  a  flow  that  has  no  permanent  cell 
boundaries.  It  consists  of  hot  rising  and  cold  sinking  transient 
bubbles  or  plumes.  At  still  higher  R  the  plumes  drift  in  one 
direction  near  the  bottom,  in  the  opposite  direction  near  the  top 
of  the  layer,  while  tilting  from  the  vertical  and  thereby 
transporting  mean  horizontal  momentum  in  the  vertical  direction. 

The  transition  when  viewed  in  the  (x,t)  representation,  is 
reminiscent  of  a  phase  transition.  When  there  are  cell  boundaries 
as  in  figure  2a,  fluid  parcels  are  forever  confined  to  circulate  in 
the  lattice-like  confines  of  the  cell  walls.  In  2b,  plumes  appear 
but  are  constrained  by  the  cell  walls  to  circulate  within  the  cell. 
In  2c  the  cell  boundaries  have  vanished  and  as  R  is  increased 
further  in  2d,  the  "melt"  begins  to  flow.  The  largest  scale  of 
motion  which  was  previously  equal  to  the  depth  d,  has  suddenly 
changed  to  the  layer  width  L.  The  appearance  of  this  new  large 
scale  has  many  implications.  Turbulent  horizontal  transport  of  a 
passive  tracer  is  completely  changed  when  this  large  scale  flow 
sets  in.  Continental  drift  and  sea  floor  spreading  which  have 
generally  been  thought  to  require  deep  convection  in  the  mantle  to 
explain  the  large  scales  observed  at  the  crust,  might  perhaps  be 
viewed  in  terms  of  a  large  scale  convective  flow  in  the  much 
shallower  but  more  plastic  asthenosphere. 
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Fig.  1.  Regime  diagram.  O  steady  cellular  flow,  •  time-dependent 
cellular  flow, [2 cells  and  transient  bubbles  (transitional),  E3 
transient  bubbles,  0m*  large  scale  flow  with  tilted  plumes,  *  heat 
flux  transitions  (see  text  and  (5)  for  further  explanation). 


Fig.  2.  (x,t)  representations  of  the  flow  0.2  cm  above  bottom 

boundary.  The  range  in  x  along  the  abscissa  is  0  to  45  cm,  t 

increases  downward,  range  in  time  is  0  to  t:  (a)  R  *  7.5  X  10^, 

P  3  0.36  X  103,  t  3  17  min.  (b)  R  3  10®,  P  3  0.86  X  103, 

t  3  22  min.  (c)  R  3  106,  P  =  7,  r  3  22  min.  (d)  R  3  2.4  X  106, 

P  =  7,  t  =  9  min. 

Fig.  3.  Sequence  of  views  of  (x,z)  slice  through  convecting  fluid, 
showing  rightward  progress  of  dye  along  bottom  and  lower-right  to 
upper-left  tilt  of  plumes. 

Fig.  4.  Heat  flux  H  vs  Rayleigh  number  R.  H  has  been 
non-dimensional ized  to  equal  NR  where  N  is  the  Nusselt  number.  The 
coordinate  scale  is  appropriate  for  the  data  shown  by  dots.  For 
the  crosses,  the  scale  should  be  decreased  by  a  factor  10.  A  least 
squares  fit  of  segments  of  data  to  a  straight  line  MR  3  aR+b  gave  a 
=  5.05  for  the  line  labelled  I,  a  3  8.40  for  II,  a  *  10.7  for  III. 
The  change  in  slope  at  a  =  1-1  X  10®  may  be  associated  with  the 
onset  of  the  large  scale  circulation  which  is  seen  at  R>2  X  106. 
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